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Electromagnetic and gravitational radiation do not propagate solely on the null cone in a generic 
curved spacetime. They develop "tails," traveling at all speeds equal to and less than unity. If 
sizeable, this off-the-nuU-cone effect could mean objects at cosmological distances, such as super- 
novae, appear dimmer than they really are. Their light curves may be distorted relative to their 
flat spacetime counterparts. These in turn could affect how we infer the properties and evolution 
of the universe or the objects it contains. Within the gravitational context, the tail effect induces 
a self-force that causes a compact object orbiting a massive black hole to deviate from an other- 
wise geodesic path. This needs to be taken into account when modeling the gravitational waves 
expected from such sources. Motivated by these considerations, we develop perturbation theory for 
solving the massless scalar, photon and graviton retarded Green's functions in perturbed spacetimes 
= ffM" + ^M"! assuming these Green's functions are known in the background spacetime g^^. 
In particular, we elaborate on the theory in perturbed Minkowski spacetime in significant detail; 
and apply our techniques to compute the retarded Green's functions in the weak field limit of the 
Kerr spacetime to first order in the black hole's mass M and angular momentum S. Our methods 
build on and generalizes work appearing in the literature on this topic to date, and lays the foun- 
dation for a thorough, first principles based, investigation of how light propagates over cosmological 
distances, within a spatially flat inhomogeneous Friedmann-Lemaitre-Robertson- Walker (FLRW) 
universe. This perturbative scheme applied to the graviton Green's function, when pushed to higher 
orders, may provide approximate analytic (or semi-analytic) results for the self-force problem in the 
weak field limits of the Schwarzschild and Kerr black hole geometries. 



I. INTRODUCTION AND MOTIVATION 

This paper is primarily concerned with how to solve 
for the retarded Green's functions of the minimally 
coupled massless scalar ip, photon Af^, and graviton 
7^2/ in spacetimes described by the perturbed metric 
gfiv = g^u + hf^i,, if the solutions are known in the back- 
ground metric g^^. One important instance is that of 
Minkowski spacetime, where these Green's functions^ 
are known explicitly. Here, we will carry out the anal- 
ysis in detail for the 4 dimensional case, and obtain 
C[ft.]-accurate solutions to the Green's functions in per- 
turbed Minkowski spacetime up to quadrature. Our 
methods are akin to the Born approximation employed 
in quantum theory, where one first obtains an integral 
equation for the Green's functions, and the 0\h^]- 
accurate answer is gotten after N iterations, followed 
by dropping a remainder term. We are not the first to 
develop perturbation theory for solving Green's func- 
tions about weakly curved spacetimes. DeWitt and 
DeWitt [1] , Kovacs and Thorne [2] , and more recently. 



Since we will be dealing exclusively with retarded Green's 
functions, we will drop the word "retarded" from henceforth. 
Despite this restriction, our methods actually apply for ad- 
vanced Green's functions too. 



Pfenning and Poisson [3] have all tackled this problem 
using various techniques which we will briefly compare 
against in the conclusions. As far as we are aware, 
however, our approach is distinct from theirs and have 
not appeared before in the gravitational physics and 
cosmology literature. 

Green's functions play crucial roles in understand- 
ing the dynamics of both classical and quantum field 
theories. The Green's function depends on the coor- 
dinates of two spacetime locations we will denote as 
X = (t,x) and x' = (t',a;'),^ and respectively identify 
as the observer and source positions. At the classical 
level, which will be the focus of this paper, it can be 
viewed as the field measured at the spacetime point x 
produced by a spacetime-point source with unit charge 
at x' . To understand this, consider some spacetime 
region V between two constant time hypersurfaces t 
and t', with t > t' . In this paper we assume that 
spacetime is an infinite (or, in the cosmological con- 



^ The spacetime coordinates in this paper will take the form 
X, x' , x" , etc. Instead of displaying the dependence on these 
coordinates explicitly, we will put primes on the indices of 
tensorial quantities to indicate which of the variables are to 
be associated with them. For example G^^/ = Gpi/[a:,x'], 
Vp'/ denotes the covariant derivative with respect to x" , g is 
the determinant of the metric at x and g' at x', etc. 
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text, semi-infinite) manifold. Let there be some field 
producing source J present in the volume V , and non- 
trivial initial conditions for the fields at i', for example, 
<p[x"^ = t'] and V° tf\x'^ = t']. Denote the scalar, pho- 
ton, and graviton Green's functions as Gx,x', G^^i and 
Gsep'tr' respectively. Then the scalar field ip evaluated 
at some point x lying on the t surface can be written 
as^ 

= / dV|5'|5G,,,.J,. (1) 
Jv 

while the photon's vector potential evaluated at x can 
be written as 



A^^ / dV|<?'|5G^,,J- 
Jv 

+ j d''-^x'\g'\^ (g^,,V°'A"^'-V°'G^,.^'^ 
In the same vein, the graviton field at x reads'* 

l^u^ I dV|g'|5G^,„-^.J' 
Jv 



(2) 



(3) 



where P'''^\,^, = (l/2)((5^,5f -f 5f 5^ - 

From equations ([!]), ^ and we see that the 
physical solution of a linear field theory can be ex- 
pressed as the sum of two integrals of the Green's func- 
tion (and its gradient), one weighted by the sources 
present in the system at hand and the other weighted 
by the initial conditions of the fields themselves. In 
particular, the d-dimensional volume integrals (with 
respect to a;') of the Green's functions, weighted by the 
field-producing Js, reaffirms the interpretation that 
the Green's function is the field of a spacetime-point 
"unit charge" because these volume integrals corre- 
sponds to calculating the field at x by superposing the 



^ These are known as the Kirchhoff representations. We refer 
the reader to the review by Poisson [4] for their derivation. In 
this paper, whenever a formula holds in arbitrary spacetime 
dimensions greater or equal to 4, we will use d to denote the 
dimensions of spacetime. Summation convention is in force: 
Greek letters run from to a! — 1 while small English alphabets 
run from 1 to d— 1 . 

* We will not be concerned with the nonlinear self-interaction 
of the gravitons in this paper. However, these nonlinear terms 
may be considered to be part of J^ji^, since gravity gravitates. 




FIG. 1: The spacetime region V is the volume contained 
within the two constant time surfaces at times t and t' , 
where t > t' . The observer is located at x = {t,x). The 
dark oval is the region defined by the intersection between 
the past light cone of x and its interior with that of the 
constant time surface at t' . On it, we allow some non-trivial 
field configuration to be present, and through the Kirchhoff 
representations in equations |T|, ([2| and (|3|, the Green's 
functions evolves it forward in time. We emphasize that 
the causal structure of the Green's function, as exhibited 
by equations Q, ([5| and ([6| means, in a generic curved 
spacetime, the observer receives fields not only from her 
past light cone (edge of the dark oval), but also its interior 
(dark oval itself) . In addition, there is some (scalar, photon 
or graviton)-producing source J which sweeps out a world 
tube, and our observer receives radiation from the portion 
of this world tube that lies on and within the interior of her 
past light cone. The picture here is to be contrasted against 
the Minkowski one, where observers only detect fields from 
their past null cone. 



field produced by all the "charges" Js present in the 
system. 

Moreover, that the Green's function yields a causal- 
ity respecting solution can be seen from the following. 
In a generic curved spacetime, if fJ^^x' (usually known 
as Sygne's world function) denotes half the square of 
the geodesic distance between x and x', a general anal- 
ysis in 4 dimensions tells us that the Green's function 
in a generic curved spacetime consists of two terms. ^ 
One of them is proportional to — i']<5[CT2:_2;'], and de- 
scribes propagation of the fields on the null cone. The 



^ See Poisson's review for the Hadamard construction of the 
Green's functions in |4|, |5| and ([6]l below. We note in pass- 
ing that, in higher than 4 dimensions, the general form of 
the Green's function will be more complicated, containing not 
only (5 and terms, but derivatives of (5-functions too. 
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other is proportional to Q[t — t']<d[(Tx.x'], and describes 
propagation on the interior of the future hght cone of 
x'. 

Gx,x' = ^ {Ux,x'6[ax,x'] + Vx,x''d[o'x.x']) , 

'ITT 

(4) 

G^.' - ^^^—^ iU^,,,S[ax,x'] + V,,,,e[ax,x']) , (5) 
47r 

_ e[t - 1'] 

'-'fj.va'P' — ~. [U fj,iya' l3'0[Ox,x' \ 

in 

])• (6) 

In Fig. pi), we illustrate the KirchhofT representations 
in 0, (If and 

Tails in curved spacetime The presence of the 
two terms in equations Q, ([s]) and ([6|, the property 
that for some fixed x', the Green's functions of mass- 
less fields are non-zero for all x both on and inside 
the future null cone for x\ teaches us an important 
difference between the propagation of electromagnetic 
and gravitational wave signals in a curved versus flat 
4 dimensional spacetime. In the latter, signals travel 
strictly on the null cone, and the radiation received at 
some location x is related to the source at retarded 
time t — t' = \x — x'\. In the former, signals travel at 
all speeds equal and less than unity.^ (We are setting 
c = 1.) This off-the-light-cone piece of massless radi- 
ation is known in the literature as the tail (or, some- 
times, wake); and is often touted as a violation of Huy- 
gens' principle in curved spacetime. As elucidated by 
DeWitt and Brehme [Jj , this implies the electrodynam- 
ics of even a single electrically charged particle depends 
on its entire past history: it exerts a force upon itself 
(a "self-force" ) , in addition to the one already present 
in flat Minkowski spacetime, because the electromag- 
netic fields it produces travels away from it but then 
scatters off the geometry of spacetime and returns to 
interact with it at some later time. 

Gravitational Dynamics This tail-induced self- 
force finds an analog in the gravitational dynamics of 
compact objects orbiting massive black holes, because 
the gravitational waves they generate scatter off the 
non-trivial background geometry and return to nudge 
their trajectories away from a geodesic one. The grav- 
itational radiation signals of such systems are believed 
to be within reach of future gravitational wave detec- 
tors, and there is currently intense theoretical work 
done to understand their dynamics. Perturbation the- 
ory in the weak field limit of Schwarzschild and Kerr 



This is barring special properties, such as the conformal sym- 
metry enjoyed by the Maxwell action, which says that light is 
blind to conformal factors of the metric: o?g^,u and g^i/ are 
equivalent in its eyes. We will shortly elaborate on this point. 



may thus provide us with approximate but concrete 
results from which we can gain physical insight from 
(and possibly serve as a check against numerical calcu- 
lations). For instance, that the tail effect is the result 
of massless fields scattering off the background geome- 
try will be manifest within the perturbative framework 
we are about to undertake; this point has already been 
noted by DeWitt and DeWitt [1] and Pfenning and 
Poisson [3]. 

Cosmology Turning our attention now to cos- 
mology, the past decades have provided us with ob- 
servational evidence that we live in a universe that is, 
at the roughest level, described by the spatially flat 
FLRW metric. In conformal coordinates, it is 

9t,v ^ a^■n^iu, '7^1^ = diag[l, -1,-1,-1], (7) 

where a tells us the relative size of the universe at var- 
ious times along its evolution. Most of our inference 
of the properties of the universe come from examining 
light emanating from objects at cosmological or astro- 
physical distances, and furthermore our interpretation 
of electromagnetic signals are based on the assump- 
tion that they travel on null geodesies. This statement 
is precisely true when the metric is Q because the 
Maxwell action that governs the dynamics of photons 
in vacuum, is insensitive to the conformal factor a?. 
Specifically, in 4 dimensional spacetime, ^Maxwell 
and S'Maxwcii[a^'7] are exactly the same object; the con- 
formal factor a? drops out. 

5'Mrxwcn[a^^] = -^^rxwciiN (8) 
= -\j d'x"7j''"''"r,-"^"F,.,„F^„p:, 

This means electromagnetic radiation in 4 dimen- 
sional spatially flat FLRW universes behaves no dif- 
ferently from how it does in 4 dimensional Minkowski 
spacetimes. In particular, it travels only along null 
geodesies. However, cosmological and astrophysical 
observations have become so sensitive that it is no 
longer sufficient to model our universe as the exactly 
smooth and homogeneous spacetime in ([?]). Rather, 
one needs to account for the metric perturbations, 

5m"^ = «^ iVfj^i' + hf^„) ■ (9) 

Because the drops out of the Maxwell action, we rec- 
ognize that a first principles theoretical investigation 
of the propagation of light over cosmological distances 
is equivalent to the same investigation in perturbed 
Minkowski spacetime. Moreover, since the geometry is 
now curved (albeit weakly so) , light traveling over cos- 
mological length scales should therefore develop tails. 
As already mentioned in the abstract, if a significant 
portion of light emitted from a supernova at cosmolog- 
ical distances leaks off the light cone, then the observer 
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on Earth may mistakenly infer that it is dimmer than 
it actually is, as some of the light has not yet arrived. 
This leakage may also modify the light curves of these 
objects at cosmological distances. To our knowledge, 
the size of the electromagnetic tail effect in cosmology 
has not been examined before. Our development of 
perturbation theory for the photon Green's function 
(and confirmation of DeWitt and DeWitt's first order 
results [T]) in perturbed Minkowski, is therefore the 
first step to a thorough, first principles based, under- 
standing of the properties of light in the cosmological 
context. This may in turn affect how we interpret cos- 
mological and astronomical observations.^ 

JWKB Now, the JWKB approximation where 
one assumes that the wavelength of the massless fields 
are extremely small relative to the characteristic length 
scales of the spacetime geometry (which, in term, usu- 
ally amounts to neglecting all geometric terms relative 
to the □ in the wave equation), is often used to justify 
that null cone propagation is the dominant channel 
of travel for massless fields in generic curved space- 
times. (See for example Misner, Thorne and Wheeler 
[B].) Here, we caution that, even in cases where the 
JWKB approximation yields exact results, it does not 
imply that light travels solely on the light cone. Such a 
counterexample is that of odd dimensional Minkowski 
spacetimes, where the momentum vector satisfies 
the exact dispersion relation rj^^^k^k^ = 0, but the 
Green's functions of massless fields develop power law 
tails: for odd d, V[x,x'] oc {{t-t'f-{x-x'f)-^'^-^^/^. 
(See Soodak and Tiersten [7] for a pedagogical discus- 
sion on tails of Green's functions in Minkowski space- 
times.) This tells us that, even for 4 dimensional flat 
spacetime, the rigorous way to prove that light travels 
on the null cone is by computing the photon Green's 
function, since it is the Green's function (via the Kirch- 
hoff representations in Q, ^ and ([3|) that deter- 
mines how physical signals propagate away from their 
sources. 

In the next section, we will review the general theory 
of Green's functions and some geometrical constructs 
related to them. Perturbation theory for Green's func- 
tions will then be delineated in the subsequent two sec- 
tions; following that, we will apply the technology to 
calculate the Green's functions in the Kerr black hole 
spacetime, up to first order in its mass and angular mo- 
mentum. We will conclude with thoughts on possible 
future investigations. 



We emphasize here that, we are not, as yet, claiming that the 
tail effect is significant in the cosmological context. Rather, 
this paper is laying down the groundwork — the computation 
of the photon Green's function - in order to investigate this 
issue from first principles. 



II. GENERAL THEORY 

This section will summarize the key technical fea- 
tures of Green's functions we will need to understand 
in the development of perturbation theory in the fol- 
lowing two sections. We refer the reader to Poisson's 
review ^ for an in-depth discussion. We first exam- 
ine the world function cr^.a;', van Vleck determinant 
Ax^x' and the parallel propagator g^^i, which are ge- 
ometrical objects needed for the formal construction 
of the Green's functions themselves. We will record 
the equations obeyed by the Green's functions, and 
then describe the coefficients of (5[cr2,_2;'] and Q[(Tx.x'] 
in Q, ([5]) and Finally we will compute the ax,x', 
Ax^x' and g^iji in Minkowski and perturbed Minkowski 
spacetimes. 

World Function The world function ax^x' de- 
fined in the introduction is half the square of the 
geodesic distance between x and x' . Assuming there 
is a unique geodesic whose worldline has coordinates 
{C"[A]|A e [0,1];C"[0] = x'", = x"}, it has the 

integral representation 

= ^^'ffM-Kl^^rdA (10) 

with i = d^/dA. 

van Vleck Determinant Closely related to ax,x' 
is the van Vleck determinant A.x,x' 

Aci[d^d^,ax.x'] ^ 

Parallel Propagator The parallel propagator 
gfj^iji is formed by contracting two sets of orthonormal 
basis tangent vector fields {e'^^jA, /i = 0, 1, 2, 3, . . . , d— 
1}, one based at x and the other at x' . (The A-index is 
raised and lowered with t^ab and the /i-index is raised 
and lowered with the metric.) 

gf,^, [x, x'] = -qABsJ^ [x\e^,^ [x'] , (12) 

with the boundary conditions that the metric be re- 
covered at coincidence x = x' , 

9t,y'[x\x']^ g^,^,[x'], gt,u'[x,x\^ g^^[x]. (13) 

The defining property of these vector fields {e^a.} ^"^^ 
hence the parallel propagator itself, is that for a fixed 
pair of X and x', the {e'^aI ^'"'^ parallel transported 
along the geodesic joining x' to x. That is, Vqe'"^ = 
and consequently 

t^o.g^..' = 0. (14) 

Green's Function Equations Next we record 
the equations defining the Green's function. For the 
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massless scalar, 



^x'Gx x' — ^xGx x' — 



(15) 



with Ux' =g'"' V^'V^' and Ux = g^^V^W^. For the 
Lorenz gauge photon (V"Aq = 0),^ 



9^.1'' 



(16) 



DeWitt and Brehme "3^ points out that the divergence 
(with respect to x) of the Lorenz gauge photon Green's 
function is the negative gradient (with respect to x') 
of the massless scalar Green's function 



V^Guiy' — l/'Gx,x' ■ 



(17) 



We will later note that our perturbative result satis- 
fies (17). Proceeding to the de Bonder gauge graviton 
(V^Tm- = with 7 = g^"^^,), 

\ {Vp} - (□ - 7^ + 2A) + 2R\^p 



{a^f3} ^ 9^^ Rap — gapR^'' jG/^^p-a 

S'^\x-x'] 



\gg'V'' 



(18) 



where we have included a non-zero cosmological con- 
stant A. The Sap-o'p' is built out of the parallel prop- 
agator g^,^'[x,x'], 

Sap-p'a' = 2 igap'gpa' + 9aa'9l3p') ■ (19) 

Green's functions are bitensors. Coordinate transfor- 
mations at X can be carried out independently from 
x' (and vice versa). Derivatives with respect to x are 
independent of that with respect to x' , so for instance, 

= d^Gap' — T^aGxp'- 

Hadamard form We now have sufficient vocabu- 
lary to describe the coefficients of 5[ax x'] and 6[(Ta;.2;'] 
in the Green's functions in equations Q, ^ and @. 
Assuming x and x' lie in a region of spacetimc where 



Our Christoflel symbol is T^, = {l/2)g°'^{d{^g,y^ - d^g^,,); 
Riemann tensor is R^^^^ = d^F'^^ + T^^F^^ - (fi ^ u); the 

Ricci tensor and scalar _R^^ = R"^^^, Tl = g^^Rpf^- Sym- 
metrization is denoted, for example, by Ts^^p^ = T^p + Tp^i- 
Antisymmetrization is denoted, for example, by Tj^^j = 
Tafi — Tfjc- Whenever we are performing an expansion in 
series of /i^i^, the metric perturbation, indices of tensors are 
to be lowered and raised with the background metric g^ii/. 



there is a unique geodesic joining them, in 4 dimen- 
sional spacetimes, the null cone pieces are built out of 
the van Vleck determinant and the parallel propaga- 
tors 



Ux,x' — \/ ^x^' 

U pu' ^ \/ ^x,x' gpu' 

U pua.' P' \f~^x^' R pv cl' p' 



(20) 
(21) 
(22) 



where 



Ppva'P' = 2 igpa'guP' + gpl3'9va' ~ 9pvga'0') ■ (23) 

The tail portions of the Green's functions satisfy 
the homogeneous equations, for example, DxVx^x' = 
Ox'Vx.x' — 0; Poisson 4] explains the appropriate non- 
trivial boundary conditions the tail function Vs must 
satisfy. Moreover, the derivation of (20 1, (21 1 and (22 1 



shows that the geometric tensors in the wave equation 
for photons and gravitons only contribute to the tail 
portion of the field propagation; while it is the differ- 
ential operator, namely □, that contributes to both 
the behavior of the null propagation and that of the 
tail piece. We will also witness this in the perturbative 
framework we are about to pursue. 

It is appropriate at this point to highlight that these 
geometrical constructs, from which the light cone piece 
of the Green's functions are built, have physical mean- 
ing for the cosmologist. For example, the world func- 
tion obeys the following equation involving the van 
Vleck determinant 



Ox<Jx 



V'cTx.K'V^ In A^^a:' =rf 



(24) 



Because S/°'ax,x' is proportional to the tangent vec- 
tor ^" at X (it points in the direction of greatest rate 
of change in geodesic distance), DxCx.x' oc Vq^" de- 
scribes the rate of change of the cross sectional area of 
the congruence of geodesies (the "expansion" ) through 
the neighborhood of x, which via (24) is related to the 
gradient of Ax^x' along the geodesies. This expansion 
scalar is related to the evolution of the angular di- 
ameter distance, which then in turn is related to the 
luminosity distance relation. (See, for example. Visser 
,8 and Flanagan et al. [S].) Along similar lines, ini- 
tially parallel null rays from an extended source be- 
come deflected due to gravitational effects (weak lens- 
ing) . Since the parallel propagator describes the paral- 
lel transport of an orthonormal reference frame along 
these trajectories, namely 



5" 



\x, X \e 



(see (12)) 



(25) 



they ought to contain physical content regarding polar- 
ization, rotation and shear of null bundles of photons. 
To sum, the light cone part of the massless scalar and 
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photon Green's function should provide an alternate 
means, from the standard ones in use by cosmologists 
today, of getting at the physics of null light traveling 
through the universe. This warrants more study. 

Before moving on to develop our perturbation the- 
ory, let us take a few moments to calculate the world 
function, van Vleck determinant and parallel propaga- 
tor up to first order in /i^i, in perturbed Minkowski 
spacetime. This will allow us to construct the null 
cone piece of the scalar, photon and graviton Green's 
function, and in turn, serve as a consistency check on 
our first Born approximation results below. ^ In fact, 
this was how Kovacs and Thorne [5j constructed the 
null cone piece of their Green's functions, by calculat- 
ing separately the van Vleck determinant and Synge's 
world function. But we shall argue that this is not nec- 
essary. The Born series scheme we have devised gives 
us a single coherent framework where all three geo- 
metric objects appearing in the null cone piece of the 
Green's function are byproducts of the computation. 
Specifically, the van Vleck determinant and the world 
function can be read off the massless scalar Green's 
function G^.x' i and the parallel propagator can be read 
off the Lorenz gauge photon Green's function G^u' ■ 

(7, A and g^^i in Minkowski The geodesic equa- 
tion in Minkowski spacetime is 



dA2 



= 0, 



(26) 



with boundary conditions i^" [0] — x'' 
The solution is 



and ^"[1] = 



i°^[\]^x'°' +X{x-x'T. (27) 
Inserting this into ( |lO| ) yields the world function 

1 



where we have defined 

A^ = {t~t\x-x'Y, 



(28) 



(29) 



which is not to be confused with the van Vleck deter- 
minant (we will always place the spacetime coordinates 
as subscripts for the latter). Since 



^^^x,x' ^^'^x,x' A^ 



(30) 



by equation (11 1, the van Vleck determinant is unity. 



In Cartesian coordinates, the parallel propagator is 



^ Some of the results here can be found in Kovacs and Thorne 
[2] and Pfenning and Poisson 3 , but we include them so that 
the discussion is self-contained. 



numerically equal, component-by-component, to the 
Minkowski metric "q^^u- 



^x,x' — l^ 



(31) 



We shall soon be making heavy use of the mass- 
less scalar Gx,x', photon G^i/', and graviton G^va'P' 
Green's functions in 4 dimensional Minkowski space- 
time, so let us record their explicit expressions here 



Gx.x' — 



Q[t - t']6[ax^x'] 
47r 

Q[t - t']5[t - t' - \x - x' 



47r|x — x'\ 

GaP' = ^apGx^x' 
G fiua' 1^' — Pfiuaf^Gx^x' 



(32) 

(33) 
(34) 



with 



Pp.,^al3 = 2 {^t^iaVnu - Vliu-nal}) ■ (35) 

The photon here obeys the Lorenz gauge ■q^^'^d^Ay — 
while the graviton the de Bonder gauge d^^h^i, — 
duri°'^hai3/2. For computational purposes, we record 
that the P has the following symmetries 



P, 



= P 



(36) 



(7, A and g^^^i in perturbed Minkowski 

To tackle these geometric entities in perturbed 
Minkowski, we start by noting that the integral in 
(10) defines a variational principle for geodesies. For 



fixed end points x' and x, and A an affine parame- 
ter, the paths which extremizes the integral in ( " 
are the geodesies. Let ^ be the geodesic in perturbed 
Minkowski spacetime joining x' to x. If we were to 
solve it perturbatively, we can try ^ = ^ -t- S^, where 
5(, can be viewed as a small displacement, and plug 
this ansatz into the integral in (10). But since the in- 



tegral defines a variational principle, that means the 
first order variation of the integrand, due to the 0[6^] 
deviation of the geodesic from the Minkowski one, is 
zero. To first order in ft,^^ , the world function can thus 
be obtained from (10) by simply setting S, — S.- 



A^A^ZW, 



where 



9(0) = 1 

^ /O 



(37) 



(38) 



(The reason for the name jj^^J will be clear later.) Now 
put (37) into (11), and employ (|30|. Then use the 



following relation, that for matrices A and B such that 
_B is a small perturbation relative to A, 

det[A + B]= det[A] (l + Tt[A-^B] + ...). (39) 



7 



(Tr denotes trace, and is the inverse of A.) We 
then deduce the square root of the van Vleck determi- 
nant is 



^/A:;^^i^l--h--h'+I^"^ (40) 



Here, h' = r^'^'^V^'' and = r]''''T^J . 

In [S] , Visser developed perturbation theory for solv- 
ing the van Vleck determinant. In particular, he 
showed that the 0[h] accurate yj'A^y is given by (his 
equation 61) 

AQA/3 rl 

v/A^«l + ^— dA(l-A)A(i?|l)„^K], (41) 

where (i?|l)Q^[^] is the linearized Ricci tensor evalu- 
ated on the unperturbed geodesic ^. 



Let us show the equivalence of (40) and (41 ). First 



we write down the explicit form of the linearized Ricci 
tensor in Cartesian coordinates. One is lead to the 
expression 



1 + ^— dA(l-A)Af9^ a^-Vv 
\d^ndfi..]\ _ \A^Af^d^d,.T!^^l (42) 



where we have employed (1 — \)d^j^ii ~ 9^' and \d^" — 
d^. In the first line, the A"' da" — d/dA. This can be 
integrated-by-parts, and the resulting (1 — 2A)9^ is 
9^ — 9^, and can be pulled out of the integral. 



A"A^ /•! 
2 X 



dA(i - \)\d^"da"hp„^„ = AP{d^ - a^')iS- 

(43) 



What remains is to demonstrate that 
A^A*^ 



dA(l - X)Xda"dp"h = -\h- hi' 



This relation can be reached by recognizing 
A°'A^da"dpi'h = d^/i/A^, followed by integrating-by- 
parts the d^/dA^. 



Equation ( 14 1 says the parallel propagator is paral- 



lel propagated along ^. If we write g^j,/ — 77^^ -I- h^^> 
and keep only the O [h] terms in the Christoffel symbol 
in Vahfj^^' = dahi_,y' - T^^l 
mately equivalent to 



ri„/iAjy', (|14|) is then approxi- 



d 

dA 



(44) 



where the derivatives are with respect to ^; for ex- 
ample, = d/d^'^. Since 6£, has to begin at 0[h], 
that means to the first order, we can replace with 
— AP. Recognizing 



dA 



APdp„h^„,„[^] 



(45) 



and recalling the boundary conditions ( 13 ) then allow 
us to integrate (44) to deduce 



1 A'' 



dit,"K"]p" K]dA. 

(46) 



As can be checked explicitly, 

dp"h^„p"[£_] = {dp +dp')K"p". 



(47) 



We may thus re- write (|46|) in terms oiiju) in (p58|). 



(48) 

To be clear, hp^ and hpi^' are the metric perturbations 
at X and x' respectively; while rip^, + hp^> is the parallel 
propagator in perturbed Minkowski spacetime. 



III. PERTURBATION THEORY 

We now describe the Born series method to solve the 
Green's functions in a formal power series in hp^, the 
metric perturbation. 

Scalar The quadratic action of the minimally cou- 
pled massless scalar field evaluated in the perturbed 
metric gp^, = gp^ + hp^ reads 

SM^\ J dV|g"|'V"Vv."^ (49) 

while the same action evaluated in the background 
metric (jp^, with denoting the covariant derivative 
with respect to it, is 

S^M cl''x"|5"l^v"VVa"^. (50) 

In 5^(5], if we replace one field with Gx,x" 1 the Green's 
function in gpu^ and the other with Gx" ,x' , the Green's 
function in g^^, upon integration- by-parts, and using 
( 15 ), we see that 



2i5'(p [ffj Gx,x" J Gx" .x' ] (51) 
- / d''x"\g"\--W""Gx,x"^a"Gx",x' ^ -Gx,x'- 
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Similarly, by replacing one of the fields in Sip [g\ with 
Gx".x' and the other with Gx^x"i one obtains 



2S'(p [(7; Gx,x" 7 Gx" .x'\ 

A\"\g"\^^"-" Gx.x"^o."Gx",x' 



(52) 



'Gx-. 



The surface terms incurred during integration-by- 
parts in (51) and (52 1 are zero because the sur- 
face integrands at hand, namely Gx.x"'^^ Gx" .x' and 

Ol — 

Gx,x"^ Gx".x', due the causal structure of the 
Green's functions, are non-zero only in the spacetime 
region defined by the intersection of the interiors of 
the past light cone of x with that of the future null 
cone of x' . As Fig. ^ informs us, this intersection is 
always a finite region of spacetime. As long as we are 
dealing with a spacetime manifold that is infinite (or 
semi- infinite) in extent, this finite region of intersec- 
tion lies deep inside the region enclosed by the surface 
at infinity, and hence does not contribute to the sur- 



face integral itself. Subtracting the equations (51| and 
(52) then hands us an integral equation for Gx,x'' 



(53) 



A\"\g"\ 



<fx" 



'S/ a"Gx.x"'^ B"Gx' 



'-9 



a" 13' 



V Gx,x"'^ a"Gx 



First Born Approximation Perturbation theory may 
now be carried out by iterating ( 53 ) as many times as 



one wishes (followed by dropping the remainder inte- 
gral terms containing Gx.x'), and expanding 



I5l 



W (1 + 2'^ 



(54) 
(55) 



to as high an order in h^j^i, as desired. (We are now rais- 
ing and lowering all indices with the background met- 
ric g^ij.) To obtain the first Born approximation, the 
©[/ij-accurate result for Gx^x'i one replaces the Gxy 
occurring within the integrals in (53) with Gx,x" and 
only need to expand the and g" ^ to first or- 

der. The result is 



Gx,x'^Gx.x'+ J dV|5"|5 (56) 

x{..^c....(l."r"^'^-.-^")a„„a...} 

with ft," = g'' hpf'r"- In perturbed Minkowski space- 
time, we set = 77^1^, employ Cartesian coordinates, 
and then use the spacetime translation symmetry re- 
flected by the Green's function Gx.x' for any d, namely 



(^fiGx.x' ^fi'Gx.: 



(57) 




V 




t' 



FIG. 2: Top Panel: The intersection of tlie interiors of 
tlie future null cone of x' and that of the past null cone of 
X always defines a finite (as opposed to infinite) region of 
spacetime. Moreover, if (and only if) x' lies on or within the 
interior of the backward light cone of x (or, equivalently, if 
and only if x Ues on or within the interior of the forward 
light cone of a;'), then there is a non-trivial intersection (in- 
dicated by the dark dashed oval) between the forward light 
cone of x' and backward light cone of x, which in 3-space 
we shall show is a prolate ellipsoid, when the background is 
Minkowski. Bottom Panel: If x' lies outside the backward 
light cone of x (or, equivalently, if x lies outside the forward 
light cone of x'), then there is no intersection between the 
forward light cone of x' and backward light cone of x. 



to pull the two derivatives out of the integral 

Gx,x' ~ Gx,x' (58) 

+ d^dp, J d''x"Gx^x" (^lh"7^^^ ~ /l""'^"^ Gx"^x' 

with h!' = r\f'^h,p"r"- This matches equation 2.27 of 
DeWitt and DeWitt [T] , if we note that their Green's 
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function is negative of ours. 

Photon Next, we turn to the photon. The 
Maxwell action in terms of electric and magnetic fields 



F^^ is 



^1 



Maxwell 



i J dV'|/|^g'^"""/'^"F^,v"F„- 



(59) 



We have already noted in the introduction, that this 
action S'MaxwoU enjoys a conformal symmetry in 4 di- 
mensions, namely, it evaluates to the same object in 
both the metric and the metric a^Qfiv] the con- 
formal factor drops out. Whenever there is such a 
conformal factor, for instance, as in the context of a 
spatially flat inhomogeneous FLRW universe described 
by the metric in ^ we will choose the Lorenz gauge 
with respect to g^i, and not a^g^j,^: 



-d, 



\9\ 



(I 



=0 



(60) 



so that the dynamics of will also be blind to a?. 
The quadratic action for the photon's vector potential 
evaluated in the metric g^^, = g^i, 



hf_,u is 



SA[g] = -\j dV|g"|5 (^^"AP"^^„Ap„ (61) 



Via steps analogous to the ones taken to obtain the 
integral equation for the scalar Green's function, re- 
placing one field with G and the other with G, we can 
write down the corresponding integral equation for the 
photon Green's function G^^' in the perturbed space- 
time g^i, = g^i, + h„ 



G 



G fJ.v' 

-x"\g 



(62) 



''/3'' A'V'X 



d''x"\g"\^(^r"^"f'"^a"G^y,Vp,, Gp„,, 
+ Gf_i\"Gp"„' 



Here and below, the barred geometric tensors such as 
Rfin are built out of g^i,; whereas the un-barred ones 
are built out of g^j, . 

First Born Approximation Like the scalar case, 
one may now pursue perturbation theory of the pho- 
ton Green's function by iterating the integral equation 
( 62 1 however many times (followed by dropping the re- 
mainder integral terms containing G^^') and perform 
the expansion in (54) and ( [55| , and of the Christoffel 
symbols 



[g] 



9] 



(63) 



)(V{„/i, 



/3}A 



to whatever order in h^^ one wishes. To 0[h], we 
merely need to replace the G^^' occurring under the 
integral sign in ( 62 1 with G^^' and develop the neces- 



sary expansion to linear order in h^^,. The additional 
complication in the photon case here, and the gravi- 
ton case below, is that one has to deal with integrals 
of the schematic form / G(r|l)VG, arising from the 
covariant differentiation of the Green's functions. The 
(r|l) is the first order in hf^i, variation of the Christoffel 
symbol, 



(64) 



For such terms, we will choose to integrate by parts, 
moving all the (single) derivatives acting on the h^^s 
in the (Fjl) onto the un-perturbed Green's functions 
G. (As already argued, there are no surface terms.) 
The ensuing manipulations require the use of equations 



(161 and (17 1. About a generic perturbed spacetime 
<?/ji^ = (jfj-v + hfj_^, we then gather that 



+ J d'x"\g"\'^ (^V„.G^v' (^lh"r"P"f'p" - h""^"f'p"^ V,.Gp..' 

+ 2^ Gi^a"h" Vq"Ga"i/' — ^Gf^cr"h°' V a"^ v'Gx" ^x' 
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1^ 



Ga"y' + a"G ^\"h'^ ^ ^ y'Gx'\x' 
"A" 



(65) 



In (651, we are again raising and lowering all in- 
dices with the background metric gp,v Here and 
below, (TV\n), {R\n)^i, and (R\n)^iyai3 are the por- 
tion of the respective geometric tensors (built out of 
5^11' = 9^ll' + ^tiu) containing precisely n powers of the 
perturbation hn,v 



When the background is Minkowski g^i, = 77^^ all 
the barred geometric tensors are identically zero. Like 
in the scalar case, we employ Cartesian coordinates 
and the spacetime translation symmetry property of 
Gx,x' in (57 1 to massage (65) into 



G fj,i/' ~ Gx^x'V^il' ~^ c^Gx^x' {^fl'u' ~t~ ^fiu) 



da; <,r]fj„ 



(66) 



M0,Gx,x" {]^h"r"^" - h''"^"^ Gx",x' 

X,x"h iyiiT.Gx'\x' "1" Gx^x" i.R\^^ ^"i/" Gx"^x' 



This matches equation 2.23 of DeWitt and DeWitt 
[T], up to a sign error, if we take into account both 
their R^p and Green's function are negative of ours. 
(Their sign error^" is the following: the two terms on 
the line right before the last line (involving the Ricci 
tensor), should both carry a negative sign each, since 
they must have come from integrating by parts the 
term —5[a\haii^iiyi6'^ W]-) As a consistency check of 
this result, one may perform a direct computation to 



show that the G^i,' in (66 1 satisfies (17) to first order 
in /i^^. 

Graviton 
Hilbert action 



Gravitation as encoded in the Einstein- 



S 



1 



EH 



167rG^ 



d'^x\g\^ {n-2K) (67) 



is a nonlinear theory. (Gn is Newton's constant and 
A is the cosmological constant.) One can insert the 



metricguy + ^/'i2T:G^^^u into the Einstein-Hilbert ac- 
tion (67) and find a resulting infinite series in 7^^. 
The quadratic piece, which will determine for us the 
Green's function of the graviton, is 



2i?.A/3M7^^7^'^ - 27^"7."i?^. + 27 l^^Rp. 
1 



7^ (7^-2A) 



(68) 



where we have chosen the dc Donder gauge V^7py = 
|V,y7, with 7 = g^'^^^v (The geometric tensors in 
(p8|, such as Ryxpfi, are built out of gfiu- 

) From mm 



and following the preceding analysis for the scalar and 
photon, we may write down the integral equation in- 
volving the graviton Green's functions 



Gsep'a' Gg^p'fjf 



dV'|g"|M Vr"Gw0"<?' 



^ k" G y." v" p' a' 
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G n" u" n' (T' 



d^ll\-ll\h ( V7 -r" K," ( ^i" -13" u" _ ^ ;^a"/9" \ vy 

2" " ) V K"'~'fi"u"p'a' 



(69) 



first iJorn Approximation Because of the number 



of terms and the plethora of mdices m ( 69 1 , the per 



turbation theory about a generic background g^j^i, and 
arbitrary dimensions d is best left for a computer al- 
gebra system to handle. We shall be content with the 
case of 4 dimensional perturbed Minkowski spacetime, 
and also set the cosmological constant to zero for now. 
To 0[h], we replace in (69) all the Gg^a'P' occurring 



under the integral sign with PseapGx.x' (see (34)) and 
expand all quantities about Minkowski spacetime up 
to first order in perturbations. Let us employ Carte- 
sian coordinates, raise and lower indices with rj^^^ , and 
integrate-by-parts the derivatives acting on ft,^^ occur- 
ring within the Christoffel symbols, 



d X G x^x" ^ ^"^a" ^" G x" ,x' 



(70) 



(a^ + a^O / d''x"Gx,x"K'-p' 

I Gx" ^x' •) 



where we have invoked (57). It helps to exploit the 



symmetries of the Riemann tensor indices {R^^ap = 
RaBp.v = -Ryfiap = -Rfiupa), thosc of P rccordcd in 
(36), and to recognize that, in d = 4 dimensions. 



(71) 



For reasons to be apparent in the next section, we shall 
re-express all the drd'^ as 



followed by using the Minkowski version of (15), 
namely 



d''Gx,x' = d'^Gx,x' ^S^'ix- x']. (72) 
We then arrive at 



Gsep'a « Gx,x' [ PStpa + | {Vpsihec + V ) + Vcysihep + h^'p') + 77pe(/lA> + ^-5 V ) + Vc^ei^Sp + hs'p')) 
~ -jVpcrhst - ^VSehp'a' 



d X { Psepadadp'Gx,x" 



Q/,"^"^ - /i""^") Gx",x' 
+ - ^^') (i^le + d[,_')Gx,x"hp'']r"Gx",x'V'7S + (5[e + d[^,)Gx'x" h„"]r 



'It" Gx" .x' 



x" .x'VpS 



+ {d[s + d[S')Gx,x"hp"]r"Gx",x'''lcre + (% + d[S')Gx,x" h^„]T.„Gx'\x'Vp( 



+ Gx,x" {PsepainlD + lJes{R\l)p"a" + Vpa{R\l)s"e" 

- \vp{5{R\'^)e"}a" - ^Va{5{R\'^)e"}p" + {R]'^) p" {S" c"}a") G x" ,x' 



(73) 
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One scattering approximation Let us examine 
(58), (661 and (73). The terms that do not involve 



any integrals can be viewed as the propagation of null 
signals, modulated by the metric perturbations mul- 
tiplying the Gx^x'- The terms involving integrals, go 

schematically as dxdx' J d*x"Gx.x"h[x"]Gx",x'- Due 
to the causal structure of the Gs, this can be inter- 
preted as a scattering process. The Gx",x' tells us our 
massless field begins at the source x' and travels along 
a null ray to x"; the h[x"] says it then scatters off the 
metric perturbations (and its derivatives) at x"; and 
the Gx,x" informs us that it then propagates along a 
null path from x" to reach the observer at a;. The full 
(scattered) signal consists of integrating over all the 
x" from which the signal can scatter off. This is the 
perturbative picture for the origin of tails of massless 
fields in weakly curved spacetime.^^ From this heuris- 
tic point of view, we can already anticipate that high 
order perturbation theory will involve more than one 
scattering events contributing to the tail effect. This 
scattering picture may also help us estimate its size 
without detailed calculations, and deserves some con- 
templation. 



IV. 5[a] AND e[a] DECOMPOSITION IN 4 
DIMENSIONAL PERTURBED MINKOWSKI 



In this section we will restrict ourselves to 4 dimen- 
sions and analyze further the first order results for 
the scalar (58), photon (|66| and graviton (73) Green's 



functions we have obtained in perturbed Minkowski 
spacetime, and show that to 0[h], concrete results for 
the Green's functions can be gotten once a single ma- 
trix of integrals (involving /iq^) can be performed. We 
will also decompose these scalar, photon and graviton 
Green's functions into their null cone and tail pieces. 
As a consistency check of our Born approximation, we 
show that their null cone pieces matches the Hadamard 



form described by equations (20), (21) and (22); this 



generalizes the analysis carried out in Pfenning and 
Poisson [3] to the case of arbitrary perturbations h^^. 

In the scalar (58), photon (66) and graviton (73) 



Green's functions results, we have to deal with deriva- 
tives (with respect to x or x') acting on the following 



We are being slightly inaccurate here, in that some of the 
dxdj.1 J d'^x"G^ j.iih[x"]Gxii terms also contribute to null 
propagation, as we will see in the next section. But we want 
to introduce thi s sc at teri ng pic ture here, because it is easier 
to see it from | |58| l, ( |66[ l and | |73[ l, written in terms of the 
Minkowski Green's function Gs, than from ( |84| l, | |86| l, and 
( |87| below, which are expressed in terms of 0|o^ j./]and the 
X-integrals in ((TSI. 



matrix integral 
1 „ 



47r' 



d X Gx^x"hci" i3"Gx" ^x' 



(74) 



with the Gx x' from ( 32 ) . Because of ( 70 ) , even the 



geometric tensor terms can be expressed as sum of 
derivatives with respect to x or x' acting on (74). For 
example. 



d''x"Gx,x"{n\l)Gx"^x' (75) 
= d^+dp+ I d^x"Gx.x"{h''"^" -v"^h")Gx".x' 



where h" is the trace of h 



and 



(76) 



In appendix ( [A| ) we show that Xafi involves the inte- 
gral of hai3 (but in Euclidean 3-space) over the surface 
generated by rotating the ellipse with foci at x and i?' 
and semi-major axis (t — t')/2, about the line joining 
X and x'. (This is the dashed oval in Fig. ([2]).) 



-^af3 [x^ X 

with 

Ial3[x,x'] 



Q[t-t']Q[ax,x'%i}[x,x'] (77) 

(78) 



1 



§2 



dn, 
— I 

47r 



-hn' 



t + t' 



X + x' 



The infinitesimal solid angle is dfl — d cos 
Cartesian components of x" are given by 



, and the 




AO 



cos 9 



(79) 



To separate the light cone versus tail pieces of the 
Green's functions, we now carry out the necessary 
derivatives on ( 77 ) as they occur in (|58|) , ( 66 ) and ( 73 ) . 



There is no need to differentiate the Q[t — f], because 
that would give us 6[t — t'] and its derivatives. Since 
this would be multiplied by either 0[(T2:,2;'] or possi- 



bly S[crx,x'], S'[<^x,x'], etc., while ax 



^AV2 < 0, 



these 6,6',. . . terms can never be non-zero when t = 
t' . Schematically, therefore, the derivatives now read 
Q[t — t']dd{Q[a]X) (where the two derivatives are both 
with respect to either x or x' or one each), which in 
turn would yield two types of terms. One is the tail 
term, proportional to Q[a\ddX and the other the null 
cone ones, proportional to either 6[a]dadX, 5[a]ddaX, 
or 6'[a\dadaX. Following that, we would impose the 
constraint Gx^x' = on the coefficients of the 6[a'] and 
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5' [a] terms. This requires that we develop a power se- 
ries in Gx.x' of Tap. Since there is at most one deriva- 
tive acting on I, however, we only need to do so up to 
linear order. (Higher order terms would automatically 
vanish once we put ct = 0.) In appendix (fAj) we find 



I, 



(80) 



The i^^p — \ Jq ha"p"[S,]dX has already been quoted 
previously in ( [38| . 

Scalar By pulling out one factor of {Att)^^ from 
one of the Gs (see (32 1), our result for the massless 
scalar Green's function in ( 58 ) can be written as 



ha"0" K]dA - 



Gx,x' ~ ^ ( S[(Tx^x' 



%dp, [\^-^^f'X-T^P^Q[ax^x']^y (81) 



with X = XP'^rjpi^. Carrying out the derivatives using 
(30) would give us, amongst other terms, the following 



oterms: 



-5' [cTx^x'] {cFx,x'X - A" A'^Iq^ 



(82) 



The first term is S[ax,x']X if we employ the identity 
— ~S[z]. The second term can be considered the 
0[h] term of S[a + A°' Al^X^p] = S[a]+S'[a]A''A^Xai3 + 



Moreover, invoking (45) and the chain rule also in- 



forms us that one of the terms multiplying (S[cr2;.j;'] is 



r(0) 



(83) 



Altogether, the Born approximation, ©[/i] -accurate 
answer, for the massless scalar Green's function may 
now be decomposed into its null cone and tail pieces 
as 



Gx,x' 



An 

+ e 



Aa A/3J(0) 



ax,x' + A"A^l(°j 



-d^'d^.x - dpd,,xp^ 



(84) 



As already advertised earlier, comparison with (37) derivatives take the form 



and (401 tells us the null cone portion of our mass- 



less scalar Green's function is indeed consistent with 
the Hadamard form in ^ and ( 20 1 . 

Photon and Graviton For the photon G^jy' (66) 
and graviton Gsep'a-' (73) Green's functions, we first 



observe that they contain respectively ry^j, and Psep'cr' 
multiplied by (84), the massless scalar Gx,x'- (Specif- 
ically, first term on the first line, and the second line 
of ([66|) for the photon; and first term on the first line. 



and third line of (73 ) for the graviton.) The light cone 
portions of these terms therefore contain the first or- 
der van Vleck determinant. For the rest of the inte- 
gral terms, we first make the observation that d^, -1-5^' 
acting on a function whose argument is the difference 
X — x' , is identically zero. The immediate corollary 



is that all the geometric terms, via (70), do not con- 



tribute to the null cone piece of the photon and gravi- 
ton Green's function because the derivatives acting on 
the & [a] leads to zero. The remaining terms containing 



^ Ap - 



(85) 



^[^x,x']^[d^ -d- ){d[^ + dy^,)X,^-,, 



where we have utilized (30) and the chain rule. Re- 
calling (46) tells us the oj^] terms on the right side 



of ( 85 1 , when added to the non-integral O [h] ones al- 
ready multiplying Gx. x' ~ i-e., the first line of (66 ) and 
first two lines of ( 73 ) - would give us the necessary 



first order parallel propagators to once again ensure 
consistency with the Hadamard form in ([5|, ( [2l| ), ([6| 
and (22 ). That is, we may now use the expressions for 
^x.x' (37), ^ Ax^x' (40) and g^^' (46) and decompose 
the photon and graviton Green's function into their 
null cone and tail pieces. To first order in /i^i,. 
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^ iJLV' ~ \ Qliv' V I^X^X'O [<Jx,x' 



(86) 



Gsep' a' 



An 



+ + dlS')Ip"]r"V'7e + + a[^02^(T"]r"»ype 



p"{S"t"}a' 



(87) 



r 



The geometric terms (i?|l)„«„^, (i?|l)ui.> and (7^|1) 



in (86) and (87) can be obtained by taking the corre- 



sponding hnearized tensors in terms of the perturba- 
tion haf3, and replacing all the hafj with and all 
derivatives da" with 3^+ ■ 



(89) 
(90) 



with da+ = da + da'- Even though these terms involv- 
ing geometric curvature are best evaluated by differen- 
tiating lap , it is necessary to record here their analogs 
to (78 1. For instance, if {R\l)ap is the linearized Ricci 
tensor, we have 

j^Q[t-t']e[a,^,,](R]l)ap 

= [ d*x"G^^^"{R\l)apG^"^x', (91) 



where 

(l^)„^[x,x'] 



(92) 



An 



(i?|l)a 



"13" 



t + t' |A| 



„ COSW, ^ 

2 ' 2 



At the first Born approximation, therefore, we see that 
a concrete expression from the perturbative solution of 



the scalar (84), photon (86), and graviton (87) can be 
obtained once the matrix integral Tap in ( 78 ) is evalu- 



ated. We also note that, suppose Zap in (78 1 has been 
evaluated; then at least when hap is time-independent 
(space-independent), there is no need to perform the 



line integral 1^^^ in ( |80[ ); rather, 1'"^^ is gotten by re- 
placing t — t' \x — x'\ {\x — x'\ t ~ t'). In such 
cases, the Born series method advocated here allows 
one to read off, as a byproduct of a single coherent 
calculation, the world function and van Vleck determi- 
nant from, respectively, the argument and coefficient 
of the (5-function in the massless scalar Green's func- 
tion; while the parallel propagator can be read off the 
coefficient of the (5-function in the Lorenz gauge photon 
Green's function. 

Gauge dependence The skeptic may wonder if 
the gauge dependence of the vector potential could 
render the tail piece of the photon Green's function 
in (861 un-physical. To that end, we note that, for 



r(0) 



fixed x' , the only pure gradient tail term in (86) is 

{l/2)dp{da — da')i°'u- Hence, the rest of the tail terms 
do not have zero curl ~ the corresponding electromag- 
netic fields are non-zero. This provides strong theoreti- 
cal evidence that the wake effect is present for photons 
propagating in perturbed Minkowski, and by confor- 
nial symmetry, in our universe too. 

Geometry and tails Let us notice that it was all 
the differentiation that took place in our work on the 
perturbative solution of the Green's functions, which 
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can be traced to the □ operator, that gave us both 
the terms in the arguments and coefficients of the S- 
functions in the scalar, photon, and graviton Green's 
function. In turn, we have identified them as various 
terms in the world function, the van Vleck determi- 
nant and the parallel propagator (in their perturba- 
tive guises). This re-affirms our assertion earlier that 
it is the differential operator □ that is solely responsi- 
ble for the behavior of massless radiation on the light 
cone. On the other hand, because of (70), at the level 



of the Born approximation, we see that the geometric 
tensors contribute only to the tail piece of the Green's 
function. 



V. SCHWARZSCHILD AND KERR 
GEOMETRIES 

As a concrete application of our formalism, in this 
section we will calculate the null cone and tail pieces of 
the Green's functions in the weak field limit of the Kerr 
geometry, to first order in the black hole's mass M and 
angular momentum S. Setting S to zero would then 
give us the first order in mass result for the weak field 
Schwarzschild geometry. These results, when pushed 
to higher orders in M and S, would provide us with 
concrete expressions for the Green's functions to inves- 
tigate the tail induced self force and more generally, 
the gravitational n-body problem, in the weak field 
limit background of astrophysical black holes. Strictly 
speaking, because S < M^, a consistent answer for the 
Green's functions would require at least a second order 
in M calculation, but since this constitutes a signifi- 
cant computational effort, we shall leave it for a future 
pursuit. 

Schwarzschild We begin with a discussion of 
the Schwarzschild case. If we choose to write 
the Schwarzschild black hole metric in (Cartesian) 
isotropic coordinates {t,x), so that there are no off 
diagonal terms, we may express 



9^1' — V^i^ ^ h^i^ 



where 



1-i 
1 + i, 



- 1 



(93) 



(94) 



M {MY {MY 
-4:^+8 — - 12 — + 
2r \2r ) \2r ) 



hij =r]^J\ 1 + 



M 
2r 



- 1 



( Ki ^{mV JmV 



0. 



(95) 



(96) 



Here r = yj Sijx'^x^ , M is the mass of the black hole, 
and we have set Newton's constant to unity, Gn — 1- 
The power series expansion of hoo and hij can be sub- 



stituted into the I-integral in ([78|. At order {M/2r f 
and beyond, the solution would of course receive con- 
tributions from more iterations and high order h terms 
from the integral equations (53), (62), and (69), and 



would likely involve two or more overlapping I-type 
integrals. Here we will focus on the first Born approx- 
imation. 

Within the one scattering approximation, the main 
technical hurdle to overcome is therefore the class of 
integrals 



1 



+ 1 r27T 



d(cos 6'")d0" 



^"^ 4tt Jo |f"[p,6l",0"| -~z\s 



1 1^+, v^+j 



(97) 



where n is a positive integer, x" has Cartesian compo- 
nents defined in (79) (so that, in particular, p = A'^ = 
t — t'), and 



X + x' 



(98) 



Because we choose our coordinate system such that 
X — x' = |A|e3, where 63 is the unit vector in the 3- 



direction, we have the following equalities (see (A4)), 
„ , , 1 



s2 — jApsin^^ cos(/)+. 



I A p sin 9^ sin ( 
lAL _ 



'+7 



s cos y+ ) = ^^63 — X, (99) 



from which we can deduce that 



s = r + r 



COSf I = 



(100) 



(The other solution (s,cos0+) — (|r— r'|, — (r+r')/|x — 
a;' I) is inadmissible because {r + r')/\x — x'\ > 1.) Here, 
r = |x| and r' = \x'\, and the azimuth angles of x and 
x' are both equal to (/)+ + tt. 

For the moment, it helps to think of x" and z as in- 
dependent vectors which we have chosen to write their 
Cartesian components in terms of ellipsoidal coordi- 
nates (p, 6" , 4>") and (s, 0+); we will also take R in 
(A5) to be simply a constant, not necessarily equal to 



The n — 1 case has been evaluated by both DeWitt 
and DeWitt [I] and Pfenning and Poisson [3] by per- 
forming a prolate ellipsoidal harmonics expansion of 
the inverse Euclidean distance \x" — z\~^. An alter- 
nate means of getting the same result is, as already 
noted by DeWitt and DeWitt, to recognize that is 
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the Columb (electric) potential of a charged perfectly 
conducting ellipsoid defined by af".^^ By definition, 
the conducting surface is an cquipotcntial one. This 
implies that the answer to has to depend on the 
s-coordinate of z only, for that would automatically be 
a constant on the ellipsoidal surface. For z lying away 
from the ellipsoidal surface, our integral must satisfy 
Poisson's equation V^i V23l(i) [s] = (with the in- 
verse metric g*^ of (A5)), which in turn is equivalent 



to the ordinary differential equation 



0= (1-^*2). 



d2l(i)[K] ^ dI(i)[K 



dK^ 



2k- 



dK 



K = s/R. (101) 



The general solution is a linear combination of 
a constant and the Legendre function Qq[s/R\ = 
(l/2)ln[((s/i?) + l)/((s/i?) - 1)]. But the asymptotic 
boundary condition implied by the integral represen- 
tation in (97) is 



lim 



Hi) 



lim 



1 



>oo 47r z 



(102) 



(When s > i?, \M\ says s/2 « |z|; s/2 essentially 
becomes the spherical radial coordinate.) The asymp- 
totic limit 



lim Qq[3/R] 



R 



(103) 



then tells us the solution for z located outside the el- 
lipsoid is 



1 



+ 1 /•2ir 



d(cos0")d(/)" 



4^7-1 7o |x"[p,0",0"]-z[s>p,^+,0+] 

s + R 



s-R 



(104) 



A conducting surface forms a Faraday cage, so for z 
lying inside the ellipsoid, the potential is position in- 
dependent and the same as that on the surface. 



1 

47r 



+1 /-aTT 



d(cos6l")d(/)" 



|f"[p,0",0"]-z[s<p,0+,0+] 



I(i)[s<p] 



R 



In 



p + R 
p- R 



(105) 



Reinstating the relationships z — —{x + x')/2, R 
|x — x'l, p = t — t' and s — r + r' , we gather 



Hi) F - a; , 

X (e[r + r'-(t-i')]ln 



(106) 



7' + r + b; — X 



+ Q[t - t' - {r + r')]\n 



r + r — \x — X 
t-t' + \x-x'\ 



t-t' -\x-x' 



For later use, let us record the following symmetrized 
spatial derivative on , keeping in mind that {di+di> ) 
acting on any function that depends on the spatial 
coordinates solely through the difference x — x' is zero: 

{d,+dS{i) = \£-x'\-^e[r + r' -{t-t')] (107) 



r + r' + \x ~ x'\ 
r + r' — \x — x'\ 



(Note that the two (5-function terms arising from dif- 
ferentiating the ld[r + r' -{t-t')] and ld[t~t' -{r + r')] 
in 1(1) cancel each other, upon setting t — t' = r + r' 
in their respective coefficients.) 

When n > 1, this conducting ellipsoid interpretation 
for the n ~ 1 case does not continue to hold; but one 
may attempt to derive a partial differential equation in 
terms of the variables (s,0+, </)+), such that some dif- 
ferential operator V acting on the kernel \x" — is 
zero. (Note that if x" and z lived in n+2 spatial dimen- 
sions, T> would be the (n -I- 2) -dimensional Laplacian, 
but implementing this scheme would involve introduc- 
ing an additional n — 1 fictitious angles and Cartesian 
components for x" and z.) The general solutions of 
this partial differential equation may either help lead 
to a physical interpretation - just as one was found for 
the ri = 1 case - or a harmonics expansion analogous 
to the one used by DeWitt and DeWitt, so that the 
resulting series can be integrated term-by-term. Be- 
cause of the cylindrical symmetry of the integral I(„), 
the final result should not depend on (f>+. We shall 
leave these pursuits for future work, and merely sum 
up the 0[A//|A|] results here. Recalling the relation- 



ship between hafi and 1^0 from (78): 



= -<5c/3MI(i) + O 



{M/\x~x'\Y 



(108) 



The Columb potential at z can be obtained by the Green's 
function type integral f d^x" ^/g!J'E[x"]/{4:n\x" — z\), where 
32 is the determinant of the induced metric on the ellipsoidal 
surface and the surface charge density S is the normal deriva- 
tive of the electric potential, S = N^di'i!, evaluated on the said 
surface. Because N^di is a unit normal, one would find that 
the combination (Px" y/g!jN^di^ is equal to the infinitesimal 
solid angle dQ in 3 spatial dimensions, up to overall constant 
factors. 



with 1(1) given by (106). 

Kerr Let us now turn our attention to a Kerr 
black hole with mass M and angular momentum S", 
with its spin axis aligned along the 3-direction.^'^ 



This 3-direction is not to be confused with the 3-direction of 
the prolate ellipsoidal coordinate system invoked during the 



evaluation of in ( |78| 
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Starting from the Kerr metric written in Boyer- 
Lindquist coordinates (see equation 33.2 of we first 
perform the fohowing transformation on tlic r coordi- 
nate 



M 
2^ 



(This coordinate transformation would yield, when 
5 = 0, the Schwarzschild metric in isotropic coordi- 
nates.) Denoting the unit vector in the 3-direction as 
63 and further define 



S = Ses, 



(110) 



to first order in both S and M, we may then write the 
Kerr metric as 



(111) 



lap, as can be inferred from (78), by replacing the 
jaf — af'l in the time argument of ha"i3" with t ^ t'; 
and replacing the t — t' in the spatial arguments of 
ha"/3" with |af— a;'|. Since the ha'ip" at hand does not 
have any time dependence, this means is given by 



(109) replacing every i — with |x — x'\ in (116). Because 



-x'\ < r- 



r' , this means the Q[t — t' — [r + r')] term 
in (1061 may be dropped and the 9[r + r' — {t — t')] 



set to unity 



14 



?(o) _ 



X In 



(117) 



World Function The world function is ax,x' 
:,x' + A^A'^I^'^'. Some calculus reveals 



2 M5. 



where 

hap [t, x] 



In a Cartesian basis, 
d 



Q/3 



S X 



dx 



S X 



dx 



S 



d d 
dx"^ ' dx^ 



1 

r 

(112) 



(113) 



1 1 \ 2(i - t')S ■ [x X x') 
r r 



M 



rr'{l + x-x') ^^^^^ 
r + r' + \x — x'\ 



r + r' — \x — x'\ 



with S ■ {xxx') = S{x^x'' 



^x'^), X = x/r, x! 



x' /r' and x ■ x' = SijX^x'^ being the Euclidean dot 
product. 

van Vleck Determinant Because the Kerr space- 
time is a vacuum solution to Einstein's equations 



The off diagonal nature^of SfJ^^ implies that the first ^ j^j^^j ^^^^^^ ^^^^ ^^^^^ ^^^^ 



order in mass Xqo and Xij for the Kerr black hole are 



identical to that of the Schwarzschild case. As for Ii 



Oil 



by referring to (74), integrating by parts the spatial 



gradient acting on r ^, and using (57) to pull the re 



suiting two derivatives out of the integral, we observe 
that it can be gotten by acting 



J^=<SX 



d_ 
dx 



d_ 
dx' 



on the n = 1 integral in (97). That is, 



(1) 



(114) 



(115) 



Altogether, to first order in mass M and angular mo- 
mentum the Kerr spacetime hands us 

T-ap = - [^M&afi + ^la^'fi^J^ 1(1), (116) 



with 1(1) given by (106). 

We will now construct the null cone portion of the 
Green's functions by computing the world function, 
van Vleck determinant, and the parallel propagator. 
From (37 1, (40 1 and (48), we recall that these objects 



may be gotten once is known. I^*^^ is related to 



gular momentum must vanish, at least away from the 
spatial origin x 7^ 0. Visser's result (41) then informs 



us that the van Vleck determinant must remain unity 
to this order, 

A,,,, « 1. (119) 

We may also confirm this by computing the van Vleck 
determinant from the world function in (118) using 
(11), or by a direct differentiation (see (40)) 



= i(%]+%'])-M(^i + i). (120) 



^■^ The following remark is in order. Because |r — r'| < — a;'| < 
r -\- r', the only way r -\- t' = t — t' = \x — x'\ can be satisfied 
simultaneously is when a null signal is sent from x' to x with 
the spatial origin (i.e. the spatial location of the black hole) 
lying on the straight line joining them (as viewed in Euclidean 
3— space), so that x ■ x' = —rr'. But we do not expect any 
signal to be able to pass through the black hole; hence, all 
terms implying such a configuration may be discarded. 
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Parallel Propagator According to (46 1, the sym- 



metric portion of the parallel propagator can be read 
off the metric perturbations, namely 



(121) 



r r 
1 



1 



At this point, it is convenient to define 



V, 



rr' (1 + X • x') 



(122) 



(One may need to recognize {r+r') — \x—x'\'^ = 2rr'(lH- 
x-x').) By a direct calculation, one may show that Vj 
is divergence- less. 



a^v,: = d' V, = 



and it also satisfies 



d d 
dx dx' ' 



0. 



(123) 



(124) 



The antisymmetric portion of the parallel propagator 
is given by 



1 



In terms of V,: , its non-zero components are 
1 



igoj' - 9jw) = {M{t - t') + A' J,) Vj 



and 



2 {9ik' - gur ) = (MA[fc ~it~ t')J[k) V,y 



(125) 



(126) 



(127) 



Tails in Kerr By recalling (107), I^i = — J7il(i) 
reads 



-\x~x'\-^e[r + r' -{t-t')] 



X Jiln 



O 



{S/\ 



x-x'\^f 



(128) 



To first order in angular momentum S, therefore, the 
Kerr spacetime does add non-trivial terms to the null 
cone portion of the Green's functions of massless fields 
in a Schwarzschild spacetime. However, the tail part 
of these Green's functions only receives additional con- 
tributions from the Kerr spacetime within the region 
r + r' >t — t' > \x — x'\ near the null cone; no contri- 
butions due to angular momentum arise deeper inside 



the null cone, t — t' > r + r' . This latter observation is 
consistent with Poisson's findings in [T^]. In fact, we 
shall find that the tail of the scalar and photon Green's 
functions are only altered by angular momentum pre- 
cisely a,t t — t' ~ r + r' , corresponding to the reflection 
of null rays off the black hole. Only the tail of the 
graviton Green's function, which is sensitive not only 
to the Ricci curvature but to Riemann as well, experi- 
ence angular momentum effects throughout the region 
\x-x'\ <t-t' <r + r'. 

Let us now proceed to compute the various pieces 
of the tail portion of the Green's functions. Equations 



84| , ( 86 ) , and ( 87 1 tell us there are only three distinct 



building blocks. Employing (116), these are as follows. 
The first is 



I(S) = 



1 



f 4:MQ[t ~ t' - {r + r')] 



+ 



2S-{xx x') 



?J/|2 



(129) 



rr' (1 -I- ir • x') 
The second is 

1 



L-S'[t-t'-{r 



(130) 



where d^± 



dfj, ± dfj,i. In terms of Vj in (122), the 



non-zero components of I^^^ are then 



?(A) 



^~Md[r + r' -{t-t')]Vj, 



(131) 



terms to zero, and 



where (124) was used to set the angular momentum 

(132) 



?(A) 



M 



2 (e[r + r'-(t-t')]V,]) 
+ J[k (^[r + r'-(t-i')]V,]). 

The third and final building blocks are the geomet- 
ric curvature terms. The non-zero components of the 
Riemann terms are 



(^|l)o.o, 



M 



%+ (e[r 



2 

M 



(<-t')]V,}), (133) 
%+Jfc](e[r + r'-(i-t')]V,),(134) 



d[k+ {Si]tQ[r 
-{i^ j) 



(135) 



where the (i O j) means one has to take the preceding 
term and swap the indices i and j. Performing the 
appropriate contractions and utilizing ( 123 ) yields the 
Ricci tensor and scalar terms 

5[t -t'- (r + r')] 



(136) 
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(7^|1) = -2M 



S[t-t' - (r + /) 



(137) 



We note that these (5-functions (the S[t — t' — {r + r')] 
and its derivative) arise from null rays scattering off 
the point mass (i.e., the black hole) at the spatial ori- 
gin. For instance, one may also arrive at ( 136 ) by 



recalling from (91) and (92) that the (-R|1)q,^ is an in- 
tegral involving the Ricci tensor over a prolate ellipsoid 
centered at (x -\- x')/2 and whose foci are at x and x'. 
Since the linearized Ricci tensor for the metric in ( 111 ) 



and (112) is 



(i?|l)„^[f] = 47r ^MSa^p + SUly (^^ ^ 5^ 
(the 47r5(^)[af] comes from —S^^didjr^^) we have 



2./ 4tt 



X + X 



with the x" in (79). This integral leads us to (136), 
if we re-express S'^^'^lx" — z\ = S[s — p]S[cos9" — 
cose+]S[(j)+ + TT - <j)"]{y/detg^j/sme")-\ using (|A6|. 

Green's Functions We may now put together the 
minimally coupled massless scalar Green's function in 
a weak field Kerr spacetime, with the geometry de- 
scribed in (111) and (112), to first order its mass M 



and angular momentum S 



G., 



Q[t-t'] 
An 



(138) 



The Lorenz gauge photon counterpart is 



47r 



'5\ 



(139) 



eK..](r;^.i(s)-fi^) + (i^)^, 



while the de Bonder gauge graviton's is 



Gs. 



ep' a' 



Q[t- t'] 
An 



We remind the reader that the van Vleck determinant 
is unity; whereas the world function <Tx.x' can be found 



in (118), the parallel propagator g^i,' components in 



(fml), and (fml), in (IT291); the components 



ofX 



A) 



pi, in (131) and (132); and the components of the 
geometric terms such as (-Rll)^^^^ in (133) through 



(1371 



As a consistency check of the building blocks 
l/i^^ and (i?|l)^j,, by employing the identities, 

5[~z\ = 5[z], 5'[~z\ = 8"[^z] ^ 5"[z\, (140) 

we have verified that the tail part of our massless scalar 



and photon Green's functions satisfy (17), or equiva- 
lently. 



= 0. 



(141) 



Causal structure Notice the I^®^ in ( 129 ) is zero 



close to the light cone, and only non-zero for late times: 
t ~ t' > r + r' . This in turns indicates the tail part of 
the massless scalar Green's function is non-zero only 
after the time elapsed t — t' equals or exceeds the time 



needed for a null ray to travel from the source at x' , 
refiect off the black hole, and reach the observer at x. 
Furthermore it is sensitive to first order spin effects 
only exactly at t — t' = r + r' . On the other hand, 
the photon Green's functions contain, in addition to 



the TpV in ([m]) and ([T32| and (i?|l)„« in ([136]); 



the graviton Green's function contain all three build- 



ing blocks, ll^u and the curvature terms in (1133 1 



through (137). TheX^^^ and curvature terms are non- 
zero only at early times |r— r'l < |af— x'l < t—t' < r+r'; 
mathematically this is because all these terms contain 
the derivative 9^+1(1) oc <d[r + r' — {t — t')]Vj . However, 
it may be worthwhile to search for a more physical ex- 
planation, for it could lead us to a deeper understand- 
ing of the tail effect. In any case, this means both the 
photon and graviton Green's functions carry non-zero 
tails throughout the entire interior of the future null 
cone of x', though their behaviors are altered abruptly 
when the time elapsed t — t' changes from t — t' < r + r' 
to t — t' > r + r' . This is because, in the former, they 
contain effects described by I^^'' and geometric cur- 
vature; while in the latter region they are, like their 
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t-t' 




FIG. 3: Time elapsed (t — t') vs. radial distance (r) view of 
the causal structure of the Green's functions in the weakly 
curved limit of Kerr spacetime. The spacetime point source 
is located at radial coordinate distance r' from the black 
hole. The dark grey area represents the early time tail 
|r — r'l < |x — x'l < t — t' < r + r'; and the light grey 
region is the late time tail t — t' > r + r' . The dashed line 
is t — t' = r + r'; the two solid black lines are t — t' = |r — r'|. 
As already noted by DeWitt and DeWitt [T] and Pfenning 
and Poisson j3], the Green's functions undergo an abrupt 
change in behavior when time elapsed transitions from t — 
t' < r + r' to t — t' > r + r' . The region t — t'>r + r' only 
receives contribution from X^^^ in (129 1, which contains 
angular momentum 5* terms only when t — t' is exactly 
equal to r + r'; while the dark grey region |r — r'| < t — t' < 



r + r' only receives contribution from X^^^ in ( 130 1 and the 



geometric curvature terms in ( 133 \ through ( 137 1. The tail 



part of the massless scalar Green's functions is only non- 
zero for t — t' > r + r', because it is entirely governed by 
Because the photon and graviton Green's functions 
depend on X^t^ and the geometric terms, they are non-zero 
within |r — r'l < t — t' < r + r' . However, their behaviors are 
altered once t — t' > r + r' , since they too become governed 
solely hyX'-^K 



scalar cousin, governed solely by I*^^^ (We illustrate 
this abrupt change in behavior of the Green's functions 
in Figure ([3]).) Finally, we observe that spin effects are 
present on the nuU cone and, in the tail, exactly at 
t — t' — r + r' , for the scalar and photon Green's func- 
tions. Only the graviton is sensitive to the full Rie- 
mann curvature of spacetime, which unlike the Ricci 
tensor and scalar, is non-zero everywhere. This is why 
the tail of the graviton Green's function contain spin 
effects within the whole region of I if— x I < t—t' < r+r'. 



VI. SUMMARY AND CONCLUDING 
THOUGHTS 

In this paper, we have developed a general Born se- 
ries expansion for solving the minimally coupled mass- 
less scalar, photon, and graviton Green's function in 
perturbed spacetimes described by the metric g^^, = 
gfiu + h^^. The key starting points are the integral 
equations for the scalar (ISSl), photon (|62|) and gravi- 



ton (69 1 cases, which were gotten from the quadratic 



portions of the actions of the respective field theories. 
From these, one performs a power series in the per- 
turbation h^i, and iterate these equations (followed by 
dropping the remainder terms) however many times 
necessary to achieve the desired accuracy. We derived 



a first order integral representation for the scalar ( 56 1 



and photon (651 Green's functions in generic back- 



grounds, and for scalar (58 1, photon (66) and gravi 



ton (73) in a Minkowski background. Furthermore, in 
(84), (86), and (87), we decomposed these perturbed 



Minkowski results into their light cone and tail pieces, 
showing their consistency with the Hadamard form. 
We reiterate that, at first order in metric perturba- 
tions, the solution of the scalar, photon and graviton 
Green's functions is reduced to the evaluation of the 
single matrix integral in ( 78 ) ; the remaining work is 



mere differentiation. Even though we have applied our 
perturbation theory only to massless scalars, photons 
and gravitons, because all we have exploited are the 
quadratic actions of the field theories involved, our 
methods should in fact apply to any field theory whose 
quadratic action is hermitian. 

As a concrete application of our formalism, we have 
calculated the Green's functions of the massless scalar 



(138), photon (139), and graviton (140) in the weak 



field limit of the Kerr black hole geometry, to first or- 
der in its mass M and angular momentum S. A subset 
of these weak field results for the Schwarzschild case 
have previously been obtained by DeWitt and DeWitt 
T, and Pfenning and Poisson '3]. Our Kerr calcula- 
tion shows that, to first order in angular momentum 
5", there will be rotation-induced corrections to these 
Schwarzschild Green's functions, only on and near the 
null cone, namely |x — x'| < t — t' < r + r' (where 
r = \x\ and r' = \x'\). Beyond that, t — t' > r + r', 
the behavior of the Green's functions changes abruptly 
and is governed solely by the mass of the black hole. 

Of the previous approaches we have studied - De- 
Witt and DeWitt [J, Kovacs and Thorne [5] and Pfen- 
ning and Poisson [5] - DeWitt and DeWitt 's seems to 
be the most general. They utilized Julian Schwinger's 
perspective that the Green's function is an operator 
in a fictitious Hilbert space, for example, G^y = 
(a:|G|a;'), from which they found its variation. How- 
ever, on the level of classical field theory, the main 
concern of this paper, our methods do not require any 
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additional structure than the quadratic action of the 
field theory at hand. Hence, we hope it is accessible 
to a wider audience. Our null cone versus tail de- 
composition was modeled after the work of Pfenning 
and Poisson [3] (except we generalized it to arbitrary 
metric perturbations), who in turn state that their 
work was based on calculations by Kovacs and Thorne 
P]. In Pfenning and Poisson's work, they wrote down 
a perturbative version of the differential equations in 
( 15 1, ( 16 ) and ( 18 1 for a weakly curved spacetime with 



only scalar perturbations and derived integral rep- 
resentations of the solutions using the flat spacetime 
Green's function Gx,x'] their methods can very likely 
be generalized to arbitrary perturbations. However, 
repeated (and un- necessary) use was made of the equa- 
tions obeyed by the gravitational potential We feel 
this obscures the fact that the solution of the Green's 
function of some field theory depends on the geometry 
but not on the underlying dynamics of the geometry 
itself. 

Cosmology We close with some thoughts on ap- 
plying our work to cosmological physics. We have al- 
ready shown that the classical theory of light in a spa- 
tially flat inhomogeneous FLRW universe is equivalent 
to that in a perturbed Minkowski spacetime. Consider 
a source of photons that turns on for a finite duration 
of time, say a gamma ray burst at redshift z = Q. We 
display in Fig. Q that not only would these photons 
sweep out a null cone of finite thickness proportional to 
the duration of the burst, but they will also fill its in- 
terior. If t is the present time, the dark oval represents 
the light that has leaked off the light cone. From our 
calculation in ( 86 1, the tail part of the Green's function 



and hence the vector potential Af,,^^^^ begins at 0[h]. 
Because the components of the stress energy tensor of 
the electromagnetic fields in an orthonormal frame T^p 
(which is what an observer can measure) is quadratic 



in the derivatives of the potential, T; 



-\dAY, 

this means deep in the interior of the null cone T^p 
itself must be quadratic in the metric perturbations 
h^v^^ In cosmology, because the metric perturbations 
h^v are believed to be sourced by quantum fluctua- 
tions of fields in the very early universe, the h^i, at a 



At the same time, we should mention that Schwinger's 1101 
initial value formulation of quantum field theory (nowadays 
known as the Schwinger-Keldysh formalism), has in fact been 
employed to tackle the post-Newtonian program in general 
relativity, itself a weak field, perturbative problem about fiat 
spacetime. See, for instance. Galley and Tiglio [11) . There is 
very likely a position-space diagrammatic calculation one can 
do to reproduce ( |58| , | |66| l and | |73| l. 
^® A consistent 0{hF\ calculation of the stress energy tensor 
that is valid everywhere, both near the light cone and deep 
within it, would therefore require the knowledge of the photon 
Green's function to 0\h?\. 




FIG. 4: At 'X', let there be a burst of photons from a 
source of finite duration. If there were no tails, these pho- 
tons would sweep out a light cone of non-zero thickness 
proportional to the duration of the event itself. Because 
of the metric perturbations hfj^u, we have shown via our 
photon Green's function calculation that light develops a 
tail in a spatially flat inhomogeneous FLRW universe. The 
dark oval represents the tail of the photon field at the 

present time t. Deep within the light cone, we argue that 
the size of the tail elTect in our universe is primarily gov- 
erned by the power spectrum of the metric perturbations, 
which is currently being probed by large scale structure 
observations. 



particular point in spacetime is a random variable and 
to obtain concrete results one would have to discuss 
the statistical average of the product of h^i, with it- 
self, i.e. {hf^^[x']hap[x"]) , the power spectrum. The 
scalar sector of this power spectrum is being probed 
by the observations currently underway of large scale 
structure in the universe, and one would have to fold 
these data into a theoretical investigation of how large 
the tail effect is in our universe. 

One way to proceed is perhaps, following Poisson 
[T^ . to start with some generic localized wave packet 
to mimic the light from a finite duration event such as 
our gamma ray burst at z = 6. We may then invoke 
the Kirchhoff representation in ([2| (with no current, 
J'^ = 0) to evolve this wave packet forward in time. At 
some later time, one can compute the ratio of energies 
in the tail piece to that still remaining on the null cone 



dct[r]i 



(142) 



This will indicate if there is a significant correction 
factor that needs to be applied to observations of ob- 
jects at cosmological distances, when inferring their 
true brightness. Because of the integrals encountered 
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in ([2) and the complicated terms in (86), however, this 
is a difficult calculation. We hope to report on this line 
of investigation in a future publication. 
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Appendix A: The Matrix I, 



a/3 



The primary objective in this section is the analysis 
of Tajj in ( 74 1 , including its behavior near the null cone 
^x,x' = 0. A similar discourse may be found in Pfen- 
ning and Poisson 0], but ours is more general because 
we performed it for arbitrary metric perturbations. 

Let us first display the integral in its most explicit 
form, using the second equality of (32): 



1 

47r 



6[t - t" - \x - x"\]5[t" - t' - \x" 



(Al) 



product Gx,x"ha" i3"Gx" ,x' , due to the causal structure 
of the Green's function, is non-zero if and only if the 
x" lie both on the future null cone of x' and on the past 
null cone of x. This can be satisfied if and only if x' lies 
on or within the past light cone of x or equivalently, if 
and only if x lies on or within the future light cone of 



If we now assume that (A3) holds, then it is the sur 



face of the ellipsoid in ( A2 1 that we need to integrate 
over, weighted by ha'/p" [t",x"]. To see this, let us em- 
ploy ellipsoidal coordinates centered at (1/2) (x + x'), 
i.e. put x" = (l/2)(f' + x)+ x'", with 



x"'[s,0,(j)] = {^Vis/^y - (i?/2)2sin6lcos0, 
\/(s/2)2 - (i?/2)2 sin 9 sin </>, 
(s/2)cos6'). (A4) 



These coordinates fix the foci to be at x and x' but 
allow the size of the ellipse to vary with s. (The 1- 
and 2-components of x'" tell us a/ — act 
as the radial coordinate in the 12-plane, and hence we 
shall require s > R. This means all volume integrals 
involve s would have limits ds.) The Euclidean 
spatial metric in 3 dimensions goes from g^j = Sij for 
Cartesian coordinates to 



diag 



(s/2)2-(i?/2)2cos2( 



S2 - i?2 

(s/2)2-(i?/2)'cos2 0, 
[is/2f - {R/2f) sin^ 
(s,0,0), s>i?. 



(A5) 



where R = |A| 
measure is 



The Euclidean volume 



We may integrate over t" immediately, so that t" 



t 



^ t' + \x" 



the constraint that, viewed in Euclidean 3-space, the 
observer at x and the emitter at x' form the foci of a 
prolate ellipsoid, with semi-major axis A°/2, defined 

by 



liateiy, so that t = i i / ^ 2 7 \ 

This in turns yields \J ^"^^[9^] = ^ ((s/2) - (i?/2) cos 9j sm0. (A6) 



Using the expressions for the components oix" in ( A4), 
we may obtain 



t-t' = \x-x" 



(A2) 



|f"-f| = J-^cos[0], 



2^1 = J + ?COS[0]. 



This implies, to get a non-zero lap-, x needs to lie in 
the future light cone of a;'. 



t-t' > \x-x'\. 



(A3) 



This means the argument in the remaining ^-function 
of the I-integrand is t — t' — s, and the -^det \x — 
= (sin 9)/2. The integral over s can be 



X —X 



For by Cauchy's inequahty, |a;— x'| < |x— x"| + |a;' — 
which means, outside the light cone t — t' < \x — x'\ < 
\x — x"\ + \x' — x"\ and no solution can be found. Fig. 
([2]) illustrates the situation at hand: we see that the 



performed immediately, and because the lower limit is 
R, it gives us Q[t—t'—R] = 6[t— t']6[tTa;^a:'] multiplying 
hai3 with s = t — t' . We are left with the angular 
integration. 



e[t- t']e[ax,x'%i: 
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1 /■ dl7, 

2 ./s2 47r 



(A7) 



i + R x + x' ^„ 



where now x'" ~ x"'[t ~ t' ,0,(j>]. 

Small fTj^.x' expansion For small ax,x', we may 
develop as a series expansion in powers of ax,x'- 
Right on the null cone t—t' = \x—x'\, and if we lie x—x' 
along the positive 3-axis, the spacetime arguments of 
hai3 take on the Cartesian components 



't + t' A° 
r^(^ + ^cos0,O,O, 2 



x^ + x'^ R 



cos 9 



(A8) 



which is equivalent to ^[cos6'] — (l/2)(x + x') + 
{l/2){x — x')cos9. This is a straight line joining x' 
to X. Let us expand about this straight line by ex- 
pressing the time component {l/2)(t + t' + Rcosd) as 



t + t' 




t-t' 



cost 



2 V V 2 y 2 

and the 3-component (l/2)(a;'^ + x'^ + A° cos (?) as 



(A9) 



x^ + x'^ 




cos 0. 



(AlO) 



Perform a change of variables in ( A7 ) cos 9 = 2A— 1 and 



Taylor expand ha/3 in powers of <7x,x' in the time and 
3-components and in powers of yJOx,x' in the remaining 
orthogonal directions. One would find it is necessary 
to expand the orthogonal directions up to second order 
to achieve a non-zero result. With ^ = x' + \(x — x'\ 



the expansion of ( A7 ) is 



1 







(All) 



Ox..' I'. /''2A-1„ , 2A-1 

dA I — dz"ha"P" - -j——jT"0"ha"i3" 







where df,, + c*!,, is the Laplacian involving only the 
directions orthogonal to x — x' . The arguments of the 
ha"p"S on the second and third lines have been sup- 
pressed; they are the same as that of the first line. 
The single derivative terms can be converted into dou- 
ble derivatives by using (45) and integrating- by-parts. 



Up to a remainder that is of 0[cr ], one can show 



cTxx' I', /2A-1„ , 2A-1 



dA 



^ dz"ha"p" —dt"ha"f}" 



d\{l-\)\{dl„-dl,)K"p.'. 



By the chain rule, 



and hence we gather, as ax 



dxf^dx'" 



{A12) 



(A13) 



0, 



1 



ha"l3"[^]dX 
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